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A kinetic theory of glasses is developed us- 
ing equilibrium theory as a foundation. Af- 
ter establishing basic criteria for glass for- 
mation and the capability of the equili- 
brium entropy theory to describe the equi- 
librium aspects of glass formation, a mini- 
mal model for the glass kinetics is pro- 
posed. Our kinetic model is based on a 
trapping description of particle motion in 
which escapes from deep wells provide the 
rate-determining steps for motion. The for- 
mula derived for the zero frequency viscos- 
ity r)(0,r) is logr)(0,r) = B-AF{T)kT 
where F is the free energy and T the tem- 
perature. Contrast this to the Vogel-Fulcher 
law log r]{0,T) = B + A/(T-T^). A notable 
feature of our description is that even 
though the location of the equilibrium sec- 
ond-order transition in temperature-pressure 
space is given by the break in the entropy 
or volume curves the viscosity and its 
derivative are continuous through the transi- 



tion. The new expression for rj (0, T) has no 
singularity at a critical temperature T^ as in 
the Vogel-Fulcher law and the behavior 
reduces to the Arrhenius form in the glass 
region. Our formula for r](0,T) is discussed 
in the context of the concepts of strong and 
fragile glasses, and the experimentally 
observed connection of specific heat to 
relaxation response in a homologous series 
of polydimethylsiloxane is explained. The 
frequency and temperature dependencies of 
the complex viscosity r]((ii,T), the diffu- 
sion coefficient D((o,T), and the dielectric 
response s((o,T) are also obtained for our 
kinetic model and fotmd to be consistent 
with stretched exponential behavior. 
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1. Introduction 



In this paper we first critically review the entropy 
theory of glasses. After defining a glass in Sec. 1.1 we 
show in Sec. 1.2 the need for an equilibrium thermody- 
namic theory of those materials that form glasses. Sec. 
1.3 gives our reasons for believing that the vanishing of 
the configurational entropy S^, or at least the entropy 
reaching a critically small value, is associated with glass 
formation. Sec. 1.4 describes briefly the many experi- 
ments that support the entropy theory of glass forma- 
tion. Sec. 1.5 offers a critique of equilibrium theories. In 
Sec. 1.6 the suggestion is made that the ^c = criterion 
can be replaced by S^ = S^o- S^o is a small critical value 
of the entropy which is dependent on the time scale of 
the experiment but is positive even for infinitely long 
time scale. Sec. 1.7 contains qualitative insights into the 
kinetics of glass formation arising from the ^c ^ 



criterion, while Sec. 1.8 makes the observation that the 
fluctuation-dissipation theorem provides quantitative in- 
sights into the connection between the equilibrium and 
kinetic properties of glasses. 

The kinetic theory is developed in Sec. 2. In Sec. 2. 1 
we pass from phase space to configuration space and 
gain an insight into the topology of configuration space. 
In Sec. 2.2 we use the principle of detailed balance to 
evaluate the transition rate constants of the master equa- 
tion describing minimal models of glass formation. In 
Sec. 2.3 using a trapping model for the phase point we 
define these minimal models and derive their associated 
(master) equations. Sec. 2.4 contains derivations of the 
zero frequency diffusion coefficient D (0, T) and com- 
plex viscosity ri{0,T), while in Sec. 2.5 frequency 
dependent D((w, r), ri{co,T) and the dielectric response 
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s{co,T) are obtained. These quantities each depend on 
the distribution of well depths W{E). This quantity 
W{E) = exp(5'c/A:) is discussed in Sec. 2.6 where viscos- 
ity is shown to be a function of free energy. Finally, 
Sec. 3.1 discusses our results while Section 3.2 offers 
some conclusions. 

1.1 Operational Definition of a Glass 

We define a glass to be a material which is an ordi- 
nary liquid at high temperatures and whose thermody- 
namic extensive quantities, volume V, and entropy S, 
fall out of equilibrium as we lower the temperature past 
some temperature T^ which depends on the rate of cool- 
ing. Above Tg the relaxation times associated with vis- 
cosity are less than the time scale of the experiment, 
while below T^ they are greater. The above definition 
describes the formation of a crystal as well as a glass so 
we augment our definition by requiring that the exten- 
sive thermodynamic quantities be continuous at Tg and 
that there be no change of spatial symmetry as we cross 
Tg. This operational definition immediately suggests a 
number of questions which must be answered if we are 
to understand glasses. 1) What are the V{T,P) and 
S{T,P) equations of state on the high temperature side 
of Tg? 2) For a given rate of cooling, why does the 
glass transition occur at one temperature, Tg, rather than 
some other temperature? 3) What are the thermody- 
namic properties well below T^ where the relaxation 
times for diffusion of molecules are so long that some 
degrees of freedom are frozen out and only oscillatory 
motions occur? Experimentally the glass is known to 
behave like an elastic solid. 4) What is the viscosity 
ri{co,T,P), where co is frequency? The first three ques- 
tions are concerned exclusively with the equilibrium 
properties of glasses. 

1.2 Necessity for an Equilibrium Theory of Those 
Materials That Form Glasses 

There are four bona fide reasons to formulate an 
equilibrium theory of glasses [1, 2]. They are: 

1) Glasses have equilibrium properties above T^ and 
well below T^. It is sensible to ask what they are. 

2) The crystal phase is not ubiquitous. This proposi- 
tion was proved in Ref. [2]. Therefore, an equilibrium 
theory is needed for the low temperature phase which 
we know is not a crystalline phase. Of course, thermo- 
dynamics is also needed to describe the low temperature 
metastable phase of those materials that can crystallize. 

3) An equilibrium theory is needed [3-5] to resolve 
Kauzmann's paradox [6,7]: An equilibrium theory 



allows us to extrapolate equilibrium quantities through 
the glass transition to see how the "negative entropy" 
and "volume less than crystal volume" catastrophes are 
avoided even when the experimental relaxation times are 
projected to be infinite. For polymer glasses the sharp 
leveling off of the experimental thermodynamic quanti- 
ties must also occur in a correct equilibrium theory. This 
either is a second-order transition or it approximates 
one. Either case allows us to calculate a T2 to which the 
Tg tends in very long time-scale experiments. 

4) An equilibrium theory is a necessary prerequisite 
for an understanding of the kinetics [7]. 

1.3 Vanishing of Configurational Entropy is the 
Thermodynamic Criterion of Glass Formation 

Once one is convinced that the equilibrium proper- 
ties of glassy materials exist there are no options. One 
simply evaluates the partition function and then the two 
equations of state V{T, P) and S{T, P). It is required, of 
course, that the important characteristics of the 
molecules be taken into account, at least within a mini- 
mal model. This minimal model (the simplest model 
which retains the essence of the problem) must have 
both intermolecular energy to allow for volume changes 
and intramolecular energy to allow for temperature de- 
pendent shape changes of the molecules. The lattice 
model of Gibbs and Di Marzio [3-5] (GD) is a minimal 
model for polymers which incorporates an intermolecu- 
lar bond energy iSh which regulates the number of empty 
lattice sites (volume) and an intramolecular stiffness 
energy As, which controls the temperature dependent 
shape changes. When this was done within the frame- 
work of the Flory-Huggins (F-H) approximation it was 
discovered that a second-order transition in the Ehren- 
fest sense was obtained and that the T{P) line separating 
the liquid state and the glassy state was given by the 
vanishing of the configurational entropy 



S,(T2,P2)=0. 



(1) 



The basic physics behind glass formation in polymers is 
as follows. At high temperatures, because of the (semi-) 
flexibility of the polymers and the large numbers of 
holes, there are many ways to pack the molecules 
together in space. At these temperatures the interfer- 
ences among the molecules are not of the kind that 
prevent the molecules from taking up their preferred 
shapes; if the internal energy associated with shape ( is 
Ei then the probability of observing shape ( is propor- 
tional to exp(-;S£',)- As we lower the temperature the 
configurational entropy approaches zero. The individual 
molecules now can no longer continue to achieve their 
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Boltzman shapes (the shapes impHed by the Boltzman 
distribution of internal energies) for as the mathematics 
show this would imply that S^ s 0, which is an impossi- 
bility. Instead the molecules are frustrated [8] by their 
neighbors from achieving their individual Boltzmann 
shape distributions and at lower temperatures (T < T2) 
characteristic of the glassy region the distribution of 
shapes of the molecules is given by the Boltzmann dis- 
tribution at 7*2 . 

Liquid Crystal Frustration 

This interpretation is strengthened by our understand- 
ing of the isotropic to nematic phase transition occur- 
ring in a system of rigid-rod molecules. At low concen- 
trations of an isotropic distribution of rigid-rods the 
entropy is large because the rigid-rods have both orien- 
tational and translational freedom. However, as the 
rigid-rod concentration increases these freedoms begin 
to disappear until at a critical concentration there is no 
longer any freedom for the rigid-rods to rotate or trans- 
late provided only that the distribution of orientations is 
random. This is the point where the configurational en- 
tropy approaches zero (there may be small pockets 
where a trapped rigid-rod can partially rotate or trans- 
late slightly). One can gain much insight into this prob- 
lem by packing pencils or soda-straws at random on a 
table-top (this is the two dimensional problem) or piling 
together rigid sticks obtained from pruning one's garden 
(this is the three-dimensional problem). It immediately 
becomes obvious that there is a critical density above 
which one can not go if the rods are to remain isotrop- 
ically distributed in space. This critical density is given 
approximately by Vx = C/x, where x is the asymmetry 
ratio of the rods and Vx is the volume fraction of rods. 
The constant C is about 4 for one lattice model [9] and 
8 for another [10]. For straight rigid-rods the system has 
a way out of the packing difficulty; the rods can align 
and do so forming the nematic phase [11]. The ordered 
phase has a larger entropy than the disordered phase 
because as the reader can readily verify by a simple 
table top experiment (partially) ordered rods gain both 
translational and rotational freedom! 

Packing of Semi-Flexible Polymers 

Semi-flexible molecules also have the option of align- 
ing. There are two cases. The first easily understood 
case is when the straight shapes are also the low energy 
shapes. In this case we form either crystals or liquid 
crystals. The second case is where the low energy shape 
is some contorted "random walk" shape. Then straight- 
ening the molecules in order to pack them in parallel 
array would raise the energy and not be preferred. 



Instead the molecules are stuck in their zero or low 
entropy contorted "random walk" unaligned state [2]. 

A Critical Entropy for Glass Formation 

The configurational entropy S^ for polymers is easily 
evaluated in the F-H approximation [3-5]. More gener- 
ally, for non-polymer as well as polymer systems S^ is 
defined as the total entropy minus the (proper extrapola- 
tion of) vibrational entropy. The volume on the T{P) 
line determined from Eq. (1) is not constant; neither is 
the number of holes in the lattice model. In fact, the 
configurational entropy can be expressed as a function, 
Scif, Mo), of the fraction of flexed bonds,/, and the num- 
ber of holes, «o- This can be seen clearly from the ex- 
pression for the partition function 2a 



Ga = E fi(f,«o)exp(-i3£(/-,«o)-i3/'y) 

/.no 



(2) 



where the volume is y = C {xn., + no),C being the volume 
of a lattice site, x the D. P., and n^ the number of polymer 
molecules. The sum is over all/,Mo such that fi(f,no) 
> 1 . Since the use of the maximum term is legitimate 
[4, 5] for this system we have S(f, «o) = klnfiif, ng). The 
condition S^(f,no) = 0, or alternatively fi(f,no) = 1, 
divides /, Mq parameter space into the large /, Mq region 
for which there are large numbers of configurations 
whose number fi(f,no) is given by exp{S^(f,no)/k) for 
each set of values /, «o and the small /, Mq region for 
which there are very few configurations because S^ = 
in this region. Both /and Mq vary along the T{P) transi- 
tion line which separates the liquid from the glass phase. 
Below T2 the values of/,no are those which obtain at T2, 
P2 when we cool at constant pressure. 

If we vary pressure below the glass temperature the 
equilibrium values of both / and «o change to those 
values appropriate to the new Tzo, ^'2,0 pair. Although 
the entropy is zero in the glassy region this only means 
that lim{SJN) = as the size of the system N ^> ^. 
There can be many allowed configurations below T^ 
consistent with this condition and this means that there 
is sufficient mobility to allow «o and /to approach their 
new equilibrium values when pressure is changed. It is 
important to realize that Mq is not a constant in the glassy 
region. Therefore, critical volume cannot be a criterion 
for glass formation. 

Since two independent equations of state (i.e., the 
PVT and the 5'VT equations) completely characterize the 
thermodynamics, within the accuracy of the lattice 
model calculation there can be no other thermodynamic 
criteria of glass formation other than the vanishing 
of the configurational entropy, Sc = 0. This important 
conclusion is supported by arguments for a relation 
between S^ and the viscosity ri{T,P) [3, 12]. The 
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physical idea for this connection is very clear. If the 
number of configurations becomes smaller and smaller 
as we approach the glass temperature from above, 
flow — which is a moving or jumping from one allowed 
configuration to another — becomes more and more dif- 
ficult and consequently the viscosity becomes larger and 
larger. This suggests that the configurational entropy 
approaching zero is the universal criterion for glass for- 
mation We now quantify the implications of the above 
statements. 

1.4 Evaluation of Tg for Polymers from the Sc = 
Condition 

If we identify the glass temperature as the point at 
which the configurational entropy equals zero then 



Sc{n,p) = 



(3) 



can be used to determine Tg. We have done this for nine 
separate classes of experiments on polymers: 

1) Tg vs molecular weight for linear polymers [1, 3]. 

2) Tg vs molecular weight for ring polymers [13, 14]. 

3) Tg vs copolymer composition [15]. 

4) Tg vs blend composition [16, 17]. 

5) Tg vs pressure [18, 19]. 

6) Tg vs cross-links in rubber [20]. 

7) Tg vs strain in rubber [20]. 

8) ACp at Tg for large molecular weight polymers [21]. 

9) Tg vs plasticizer (diluent) content [22, 23]. 

In all cases we obtain reasonable fits to the experi- 
mental data. There are several interesting aspects to 
these comparisons. First, there are essentially no 
parameter fits to experiment since the model parameters 
are determined by other independent measurements. In 
item 1) of the above list we fit to the glass temperature 
at infinite molecular weight in order to determine the 
stiffness energy As (one parameter). In 5) we need to 
assume how the volume of a lattice site varies with 
pressure (one parameter), the remaining theoretical pre- 
dictions involve no parameter fits to experiment. 

Each class of experiment illustrates a feature of poly- 
mer glasses. Item 9) illustrates the colligative-like prop- 
erties of glasses. The initial glass temperature depres- 
sion by low molecular weight diluent is predicted [23] 
to obey the equation 



jdTN/dm = -3Tg 



(4) 



where m is the total mole fraction of diluent expressed 
in terms of mole fraction of monomers, and y is the 
number of flexible bonds per monomer. One notices the 
universal character of the prediction. Item 5) predicts 
that Tg vs pressure curves have horizontal asymptotes at 
high pressure. On the other hand, the free volume 
theory which assumes that the glass transition occurs 
when the hole fraction reaches a critically small value 
(usually 0.025) predicts a vertical asymptote. 

In 8) the specific heat change at Tg for a large molec- 
ular weight polymers is given to within 10 % by [21] 

ACp = Rfil-f)iAs/kTgf RTgAai4-TgAa/0.06) 

+ 0.05TgAaCj,iTf), (5) 

where R is the universal gas constant,/is the fraction of 
flexed bonds at Tg, Aa is the change in the thermal 
expansion coefficient as we pass through the glass tran- 
sition, and Cp {Tg~) is the specific heat just below Tg. 
Notice that this is a no parameter prediction since Tg, 
Aa, Cj,{Tg~) are known from experiment, As/kTg is 
determined from the condition that 5'c(T'g) = 0, and /is 
a known function of only As/kTg (f= 2exp{- As/kTg)/ 
[l+2exp(-Ae/A:T'g)]. In 2) the glass temperature is pre- 
dictcf^ to rise as we lower molecular weight of rings in 
accordance with experiment. This is purely an entropy 
effect [13, 14] arising from the observation that a ring 
of molecular weight x has more entropy than two rings 
each of molecular weight x/2. Thus a bulk system of the 
larger rings, since it has the larger entropy, must be 
cooled further to reach the Sc = condition which 
defines Tg. It should be noted that the fits of theory to 
experiment have all been made with the original Gibbs- 
Di Marzio theory [4, 5]. We have not needed to adjust 
the theory to account for new experimental data. 

Finally, we should remark that a perfect fit to exper- 
iment would require that a) the F-H calculation is 
perfect. It is not, because the statistics are approximate 
and because the molecules are modeled imperfectly; b) 
that the experimental data is excellent, including the use 
of well characterized polymer material; c) that kinetics 
have no sensible effect on the comparison with experi- 
ment. We would argue that since kinetics are important, 
perfect accord with experiments would be proving too 
much. We are predicting the underlying transition tem- 
perature T2, and the relation between T2 and the experi- 
mental Tg requires further elucidation. We should note 
that our theory predicts T2 and not Tg. Since in our 
equations T2 appears only in the dimensionalless forms 
As/kT2 and EJkT2, if the predictions are correct for T2 
and if T/Tg depends only on the rate of cooling, then the 
predictions for Tg will also be correct. Our good fits to 
experiments suggest that Tz/Tg is a constant or T2 and 
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Tg are not very different; or some combination of the 
two. A question we have not examined is "If the criterion 
for glass formation is ^c -^ S^,o how well does it predict 
glass temperatures?" It may suffice for ^c.o to be small 
(see below). Mention should also be made of the 
attempts to predict the glass temperature of a material 
by simply noting the chemical structure. Figure 10 of 
Ref. [17] and Fig. 6 of Ref. [24] are remarkable and 
suggest that further progress can be made. In both of 
these predictions an entropy criterion is used. 

1.5 Critique of the Correct Equilibrium Theory of 
Glasses 

An equilibrium theory must satisfy the following 
criteria: 

1) Accurate predictions of thermodynamic quanti- 
ties without multiplication of parameters. 

2) It must explain the ubiquitous nature of glass 
formation. 

3) It must explain why glasses fall out of equi- 
librium as the glass temperature is approached from 
above. 

4) All predictions must be correct. Since the lattice 
statistics used for glasses are applicable without change 
to rigid rod molecules, if the theory is applied to liquid 
crystals the predictions for this class of materials must 
also be in accord with experiment. 

5) It must provide a foundation for kinetic theory. 

We believe we have done reasonably well with regard 
to criterion 1) as the previous section indicates. 

Criterion 2) may be met by first defining the config- 
urational entropy for all materials as the total entropy 
minus the extrapolation of the vibrational entropy. Any 
method of evaluating the partition function from first 
principles which gives the proper equilibrium behavior 
above T^ is viable. One would then identify the glass 
transition as the place where ^c becomes smaller than 
some critical value as we cool the system. The following 
systems need to be examined for their glassy behav- 
ior: (a) Polymer glasses, (b) low molecular weight 
glasses, (c) The classic inorganic glasses, (d) liquid 
crystals, (e) systems composed of plate-like molecules, 
(f) spin glasses, (g) plastic crystals, (h) metallic glasses, 
and (i) gels and thixotropic materials. A common fea- 
ture of these diverse materials is that they each show 
frustration-the molecules, or spins, are each prevented 
from achieving their preferred low energy shape by the 
interferences of their neighbors. See below. 

Under 3) we must be careful not to equate falling out 



of equilibrium with loss of ergodic behavior. There is a 
sense in which a system is never ergodic, even at high 
temperatures. To see this for the case of polymers 
consider a polymer of A^ monomers which we model as 
a self-avoiding walk (Sj'W). An estimate of the number 
of configurations of one polymer molecule, assuming a 
cubic lattice is given [25] by 4M"N"^ » 4" = 10."'^" 
For A^ = 1000 which is a small molecule for polymers, 
the total number of configurations that can be sampled 
during the lifetime of the universe which is about 
10'" years is 10^^ moves/s X 1000 monomers X 3.6 X 
lO's/yr X 10'" yr = 3.6 X 10^^. This number is so much 
smaller than 10^"" that we see immediately that no sys- 
tem is ever ergodic. Obviously, effective ergodic behav- 
ior over some time interval is the relevant concept. By 
falling out of equilibrium we mean nothing more than 
that there are certain correlated motions of the 
molecules that occur with less frequency as we cool the 
system. At the glass temperature and below they are so 
rare as to be not measureable. In general one expects 
that the glass temperature depends on the particular 
correllated motion being used to monitor it as well as 
the rate of cooling. 

Under 4) above we have the happy circumstance that 
the same F-H lattice model that was used for glasses 
also predicts the formation of liquid crystals. As On- 
sager originally observed [11], the nematic phase of 
liquid crystals occurs because of the increased difficulty 
of packing rigid rod molecules together in space as we 
increase their concentration. Thus, the isotropic to 
nematic transition in liquid crystals occurs because it is 
entropy driven — configurational entropy driven. The 
nematic liquid crystal phase occurs for the same reason 
as glasses and the correctness of the F-H calculations for 
liquid crystals argues for their correctness for glasses, 
and conversely. The transition from random order to 
parallel alignment for a system of plate-like molecules is 
also entropy driven [9]. Although the decrease in con- 
figurational entropy drives the transition in all three 
cases the results are somewhat different. Rods and 
plates have a way out of the packing difficulty; they can 
align, thereby increasing the configurational entropy. 
Rods lying in parallel with some freedom about the 
director have a higher configurational entropy than a 
random packing of rods that is up against its dense 
packing limit. The molecules forming a glassy material 
may not have this option. To see this, suppose the lowest 
energy shape is chosen to be such that the molecules, if 
they each have this shape (and if we specify that the 
packing leaves no lattice sites unoccupied), cannot pack 
in regular array on a lattice; the majority of polymer 
shapes are of this type [2]. Then alignment at low tem- 
peratures is not favored and the material is stuck in its 
glassy phase. 
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Other Entropy Theories 

It is important to improve on theoretical predictions 
of the equilibrium properties of glass forming materials. 
One can not expect that the Gibbs-Di Marzio theory 
which is an elaboration of the Flory-Huggins lattice 
model is the final word. Improved equations of state 
would permit more critical tests of the entropy hypothe- 
sis to be made. An improved theory should derive the 
P-V-T and 5- V-T equations of state to equal accuracy. 
A theory that gives a poor 5 -V-T equation of state is sure 
to give undue stress to imagined implications of the 
P-V-T equation of state. The theory must allow the 
molecules to have shape dependent energies, since these 
are undoubtedly very important to glassification in poly- 
mers. We stress that an improved theory may not show 
an actual underlying second-order transition as ours 
does (there may be a rounding), but it should approxi- 
mate one. 

Two theories that include the effects of stiffness en- 
ergy are those of Gujrati and Goldstein [26] and of 
Milchev [27]. We accept that Gujrati has calculated a 
rigorous lower bound to the entropy for a two dimen- 
sional square lattice. This means that we must modify 
our criterion of glassification, viz. ^c — > 0, to something 
else. 

We do not accept the Milchev criticism [27] because 
his formula does not show the phenomenon of frustra- 
tion which we take to be an essential feature of glassifi- 
cation and, for rigid rods, an essential feature driving the 
isotropic phase towards the nematic phase. Specifically, 
in the Milchev theory individual polymer chains are 
never prevented from achieving their Boltzmann 
distribution of shapes which are given in the simple 
nearest neighbor model by / = (z-2)exp(-/A:r)/[l + 
(z - 2) exp(- AeAT)] where z is the coordination number 
of the lattice. In the Gibbs-Di Marzio model this distri- 
bution is realized above T2, but at lower temperatures 
each chain is frustrated by its neighbors from achieving 
its Boltzmann distribution of shapes. Instead, the distri- 
bution that existed at T2, P2 persists as we lower the 
temperature at constant P2. The number of holes also 
remains constant below T2 while in the Milchev theory 
it continues to decrease. The fact that experimentally the 
volume versus temperature curve for a glass parallels the 
volume versus temperature curve for a crystal supports 
the view that the number of holes is constant below Tg. 
It must be mentioned however that some computer cal- 
culations exist that support the Milchev formula [28]. 

1.6 Modification of the 8^ = Criterion to S^ = Sc,o 

One reason for the configurational entropy to be 
somewhat greater than zero at the glass transition has to 



do with the concept of "percolation of frustration" as a 
criterion of glass formation. As an entree to this problem 
we express the configurational entropy as a function 
Sc(f, no) of the two order parameters / (the fraction of 
flexed bonds) and Mq (the number of empty lattice sites). 
The equation 



S,(f,no) = k\nn(f,no) = 



(6) 



divides/, Mo parameter space into two regions, the large 
/, Mo region being the liquid region and the line defined 
by Eq. 6 gives /,Mo values appropriate to the glass. 
Because there are spatio-temporal fluctuations in / and 
Mo, if we are in the liquid region just above the glass 
region there will be "clusters" of polymer for which the 
f,no values are appropriate to the glassy state, and clus- 
ters for which the /,% values are appropriate to the 
liquid state. As we lower the temperature these glass- 
like clusters grow until they span the space or percolate. 
However, as is characteristic for percolation [29] there 
will be pockets of liquid-like clusters (regions of mate- 
rial for which the f,no values are appropriate to the 
liquid phase). The glass temperature would be defined 
as the highest temperature for which there is percolation 
of the glass-like structures. Because of the existence of 
the liquid-like pockets this Tg would correspond to a 
configurational entropy somewhat greater than zero. 
This percolation view of glasses receives support from 
experiments which show anomalously high mobility as 
the temperature is decreased through the glass transition 
[30]. The unexpectedly high mobility seems to arise 
from pockets of fluid dispersed in a glassy matrix. 

The percolation argument can be quantified by allow- 
ing / and Mo to vary in space. Consider the following 
diagram (see also Fig. 1). 




123456789 10 

Here the numbers indicate positions in space (cells) 
and the dashes above a given number denote values of 
/(r),Mo(r). The line connecting the 10 places is one 
particular enumeration of /,Mo in space. Each cell 
contains n lattice sites and there are N/n cells where A^ 
is the total number of sites on the lattice. To obtain the 
total partition function we take the product of the 
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partition functions for each cell. Thus, 

Q = nq(r) (7) 

where the product is over space, and 

q(r) = Xw(f(r), no)exp(-l3u(f(r),no)) (8) 

where the summation is over all permitted values of/, no 
at the cell labeled by r values on the high side of the line 
S (f, «o) = 0) in/, Ho parameter space. It is evident that all 
possible/(r), no(r) values are thereby accommodated. It 
is easy to see that the value of entropy calculated from 
this procedure is larger than that calculated from Eq. (2) 
for the simple reason that Q includes a sum over many 
paths while Qa does not include such a sum. This argu- 
ment suggests that we should replace the Sc ^ crite- 
rion for the occurrence of the glass transition by the less 
stringent criterion 



So(T2,P2)^S,,0 



(9) 



where Sc.o is some critical value of the configurational 
entropy. This is in accord both with the ideas of perco- 
lation [29] and with the experimental observation that 
small pockets of polymer within the glassy region can 




Fig. 1. In the lattice model version of the entropy theory of glasses 
there are two global order parameters:/, the fraction of bonds flexed 
and Vo, the fraction of empty lattice sites. The glass transition occurs 
when these quantities which decrease with decreasing temperature 
and increasing pressure have values which make the configurational 
entropy Sc(f, Vo) equal to zero. However, since /and Vo have spatial 
and temporal fluctuations we know that as we approach the glass from 
above there will form pockets of material for which the order parame- 
ters are appropriate to the glass imbedded in a sea of liquid with 
regions of/ and Vo appropriate to the liquid state. Percolation theory 
tells us that when these pockets connect up into an infinite cluster 
there will remain pockets of liquid. If we define the thermodynamic 
glass transition as the percolation point then the configurational 
entropy will be greater than zero at the transition temperature. 



show kinetic behavior that is not simply vibrational 
behavior (crankshaft motion [31], for example), ^co 
would be that value of entropy for which the glassy 
regions first percolate as the temperature is lowered. 

1.7 Qualitative Insights Into tlie Kinetics of Glass 
Formation Arising From the Sc^> Criterion 

We wish to determine how the kinetics relates to the 
configurational entropy Sc as we cool our system. Now, 
how one speaks about entropy depends on the kind of 
ensemble one is working with. We have used [4, 5] the 
Canonical Ensemble for which SJ^Vk = 2/ In/ but 
because the system size is large we can also write ScCr)/ 
k = lnW(T) where W(T) is the number of configurations 
whose energy is the average energy E(T) determined 
from the Canonical ensemble. This enables us to speak 
in terms of the microcanonical ensemble. 

Thus, as we lower the temperature the number of 
configurations W decreases so that they are farther apart 
in phase space [that part of phase space for which the 
total energy is E(T)]. The process of diffusion as well 
as the process of flow can be viewed as a jumping out 
of a deep well, a subsequent wandering about the phase 
space between deep wells and a dropping into a deep 
well different from that from which it had exited. The 
process then repeats itself. Obviously this process be- 
comes more infrequent at lower temperatures resulting 
in increased viscosity and decreased diffusion. There 
are several reasons for this. First, the deeper the well the 
longer the time to escape it — wells are effectively 
deeper as we lower temperature; second, the further 
apart the wells are the more time it takes for a phase 
point to wander from one well to another; third, the 
further apart the wells the larger the probability that the 
wandering phase point will fall back into the well it has 
just escaped resulting in no net flow. See Appendix A 
for a discussion of this effect. The picture we are using 
is a variant of the trapping model with the difference 
that instead of an atom or an electron being trapped we 
are trapping the phase point (or configuration point 
[32]). In Sec. 2 we shall quantify these ideas. 

1.8 An Insight From the Fluctuation Dissipation 
Theorem 

Whenever there is a thermodynamic phase transition 
the fluctuation-dissipation (F-D) theorem [33] suggests 
that dissipative quantities have the same discontinuities 
as the underlying thermodynamic phase transition: A 
simple example of a F-D theorem is the Green-Kubo 
[34] relation 



/ 



D = (l/3) (v(0)-v(O)d? 



(10a) 
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which relates the diffusion coefficient D to the autocor- 
relation function of the particle velocity v . More gener- 
ally the F-D theorem relates x (defined as the response 
of a material at {r,t) arising from an impulsive force at 
(0,0)) to the correlation in fluctuation at these two 
space-time points [25]. Since the fluctuations of a sys- 
tem at equilibrium show a discontinuity of the same 
character as the thermodynamic extensive variables, so 
do also the dissipative quantities. Thus, for a system 
undergoing a first-order liquid to crystal transition the 
viscosity 7](w,T,P) will show a discontinuity as a 
function of T,P since the volume and entropy do. 
Similarly, for a system undergoing a second-order tran- 
sition we can expect that the viscosity will show discon- 
tinuities in slope since the volume and entropy do (see 
later in this paper however). There are many examples in 
the literature of dissipative quantities such as viscosity, 
diffusion coefficient, electrical conductivity, particle 
conductivity and thermal conductivity which show 
breaks as a function of temperature as we pass through 
the glass transition. However, it is also true that a 
genuine falling out of equilibrium will also cause the 
same kind of behavior. It is uncertain how one distin- 
guishes between the two effects. Movement of the tran- 
sition point as a function of the time-scale of the exper- 
iment seems not to be a distinguishing characteristic 
since this happens also for systems known to have 
genuine first-order transitions-supercooling being an 
obvious example. 

More generally the frequncy dependent diffusion 
coefficient is given by 



-i 



D((o, T) = exp(+(wO (v(0) v(0) dt . (10b) 

•'0 

2. Kinetic Tlieory of Glasses 

2.1 A Remark on the Topology of Phase Space 

The potential energy surface of a liquid 
E'(. ■ qj. ■) appears in the partition function Q 

Q=t exp(-(K({..qj,pj..}) 

■'[■■ij.Pj--] 

+ E\. .qj. .))lkT) Ildqj Ildp, 

= A'^\ exp(-£({ . .qj.. })lkT) ndq, , (11) 

where q-^, pj are the generalized position and momentum 
coordinates of the N particles, K is the kinetic energy, E' 
is the potential energy and A the thermal wavelength. 



Since the kinetic energy is quadratic in /?, the integration 
over/7, is straightfoward. In polymers, even if £' is pair- 
wise additive, E is not [35] because the coefficients of 
the quadratic terms in K are in general dependent on qi . 
The simplification of Eq. (11) allows us to work exclu- 
sively in configuration space. This is generally repre- 
sented as a multi-well potential energy surface. As one 
approaches Tg from above the wells effectively become 
very deep because of the MkT term. One then talks 
about flow as a motion from one deep well to another 
deep well via the higher energy continuum. 

Here, however, we wish to emphasize a different 
aspect of the phase space topology. Consider the config- 
uration space of N identical noninteracting hard point 
particles on a line of length L . The partition function is 
given by 



Q = L^/Nl 



(12) 



The volume of configuration space for this system is 
given by L". Now consider the case where the particles 
each have a diameter d . The partition function is 



Qi = {L-NdflN\ 



(13) 



The ratio of the two phase space volumes is given by 
QdlQ = (1 -Nd/Lf < expi-N'-d/L) 



= e.xp(-4>N) , 



(14) 



where (p, being the volume fraction occupied by the 
particles, is on the order of 1 . Since N is on the order of 
Avogadro's number we see that the fraction of the 
volume of phase space occupied by the extended parti- 
cles is infinitesimally small relative to the unconstrained 
particles. Based on this picture a point in configuration 
space wanders on the finest of gossamer threads [36] 
which pervade the A^-dimensional hypercube of phase 
space as a fine network whose total volume is an 
infinitesimal fraction of L". The application to glasses is 
in the observation that as we lower T the effective value 
of d increases, resulting in even fewer and finer 
gossamer threads for the phase point to travel on. Thus, 
not only are the number of paths (threads) between two 
phase points fewer as we decrease T but also as one 
traverses a given thread the potential energy minima are 
effectively deeper and the barriers effectively higher. 

The above discussion serves to show how important it 
is to know how the deep wells are connected to other 
wells. In order to solve this problem we need to con- 
struct a model of the topology of configuration space 
and to calculate the transition rates for jumping from 
well to well within this model. 
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2.2 Detailed Balance Makes A Significant State- 
ment Concerning the Kinetics of Glasses 

Boltzmann's law gives sx.p(- EJkT) as the fraction of 
time that a system spends in state ( but it does not say 
how often the system jumps from state ( to state j. To 
determine this we use the principle of detailed balance 
in the form 

NiUij = NjUji , aijiaji = NjINi 

= expi-[Ej - EiVkT) , (15) 

where ay is the rate of jumping from state ( to j and A^, 
is the fraction of time a system spends in state ( . In using 
Eq. (15) one must first decide how the energy is appor- 
tioned into forward and backward transitions. For deep 
wells it is sensible to assume that all of the barrier is in 
preventing the phase point from jumping out of the well. 
It does this at a rate given by 1/t where t is the average 
time to exit the well. If we also recognize that the prob- 
ability of jumping out of the well is exponential in time 
[37] we have 



P{t,T)= T~'cxp(-t/T) and 
T-' =b,j cxpi- [Ej-E,]/kT), 



horizontal bonds of Fig. 2 is compounded from those of 
Fig. 3. 

— o p- 



K, 



A, 



<v_ 



M, 



(16) 



Fig. 2. Our minimal model for describing the kinetics of glasses. The 
points are points in configuration space and the connecting lines 
represent allowed transitions between points. The horizontal lines 
with rates 05 for traveling to the right and /3, + 1 for traveling to the left 
represent travel of the configuration point among the "configurational 
sea" of shallow wells. The vertical lines connect the "configurational 
sea" to the deep wells, the length of the vertical line being propor- 
tional to the potential energy depth of the well. The rate of escape 
from the deep wells is Aj and rate of capture is b . When the configu- 
ration point is in a deep well there is no motion; motion occurs only 
when the configuration point is cruising the configurational sea of 
shallow wells. This trapping model allows us to infer an important 
contribution to the complex viscosity Tq*(co, T), the diffusion coeffi- 
cient D(ft), T) and dielectric response e(ft), T). 



where P is the probability density of exiting the well at 
time t. It is imagined that once the phase point has 
escaped the well it wanders around in the configura- 
tional sea of the high energy region of phase space until 
it falls into a low lying well, starting the flow process all 
over again. This configurational sea consists of many 
shallow energy wells, so it is expected that jumping out 
of the deep wells are the rate determining steps. 

2.3 Diagrams for Our Minimal Models and Their 
Associated Equations 

The diagram for our primary minimal model is dis- 
played in Fig. 2. This diagram is a contraction of a vastly 
more complicated diagram but we believe it retains the 
essential features of glassy behavior. The points on the 
upper line represent the multitude of shallow wells 
while the horizontal lines connecting these points repre- 
sent the transition rates between these wells. This set of 
horizontal lines and points represent the vastly more 
complicated diagram of Fig. 3. At high temperatures 
this "configurational sea" of shallow wells is where all 
the action is; The configuration point jumps rapidly 
from well to well. The occupation number Nj for well j 
in Fig. 2 is really the sum of the occupation numbers 
vertically above it in Fig. 3, and the transition rate for the 



,'f\ ^i< ^<\ ''\- ^,' /\ ■/ 

: ■^.^ V'' I v 1 V ' V \/ » / 



^. I ^^ 




Fig. 3. The set of horizontal lines and their connecting points in Fig. 
2 really represent the vastly more complicated diagram of Fig. 3. The 
occupation probability Nj of Fig. 2 is really the sum SW,-,- of Fig. 3 and 
the a and /3 of Fig. 2 are compounded from the rate constants of 
Fig. 3. The net result is that the a and /3 are much larger than h and 
A in Fig. 2. 



The lower points represent the deep wells. Our view 
of what happens is as follows. At low temperatures the 
configuration point is in one of the lower wells. After a 
long period of time it jumps out and wanders about the 
configurational sea of upper wells until it falls into a 
low lying well. It then stays in this well for another long 
period of time until it jumps out repeating the process, 
and so on. The situation at high temperatures as 
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described in the preceding paragraph is very different. 
There are so few deep wells relative to the number of 
upper wells that they are unimportant; all the motion is 
jumping among the upper wells. The rate constants for 
jumping out of these lower wells are much much smaller 
than that for jumping back down into the well and than 
those for traveling horizontally. By adjusting the ratio of 
the rate constant for falling back into the deep well to 
that for traveling horizontally we can control the 
accessibility that the configuration point in the configu- 
rational sea has for the deep wells. 

The length of the vertical line connecting the deep 
well to the upper well(s) is proportional to the well 
depth. These vertical lines represent many possible 
paths in configuration space leading to the deep well. In 
Fig. 4 we have listed some of the possibilities. Figs. 4c, 
4d, 4e can each be shown to be equivalent to Fig. 4b. To 
see this, one writes down by the methods of Ref. [38] 
the set of equations corresponding to a given figure and 
then one shows that they can be transformed to the set 
of equations describing Fig. 4b. The rate constants in the 
transformed set of equations are such that the occupa- 
tion probabilities at each level are the same as those in 
the untransformed figure. 



<5 
a 



U 




However, Fig. 4f has a different structure entirely. In 
a descent of the configuration point from the configura- 
tional sea into this structure it can get hung up in a 
branch so that it may take a long time for it to reach 
equilibrium. The other figures all equilibrate rather 
quickly. 

The results of this paper will allow us to conclude 
that although Fig. 3 is rather simple it does catch the 
essential features of glassification. 

The Master Equations describing the minimal model 
of Fig. 2 are given by the simple set of equations 

dNJdt = - (uiNi - P2N2) - btN, + A,Mi 



dN,/dt = (a,_,Ar,_, - m) - {a,Nj - A^iA^;..) 



-b,N,+AM 



dMj/dt = + bjNj 



-AM 



(17a) 



(17b) 



where the Greek symbol rate constants denote stepping 
to the right (a) or left (;S) and the Roman symbol rate 
constants denote stepping down (b) or up (A). 

2.3.1 Going from Phase Space to Configuration 
Space to Real Space We have already shown in Sec. 
2. 1 that one can integrate over all the momentum vari- 
ables of phase space so that we deal only with position 
variables (configuration space). We would like to go 
further and deal with the smallest number of position 
variables possible. We begin by supposing that there are 
two separate noninteracting regions of space each with 
their own master equations 



Fig. 4. The vertical lines in Fig. 2 represent many possible paths in 
configuration space leading to the deep wells. In Fig. 4 we have listed 
five possible paths to deep wells, or equivalently ways to decorate 
each of the vertical lines of Fig. 2. It can be shown that diagrams b, 
c, d, and e are equivalent to a. Thus, the diagram of Fig. 2 really 
represents vastly more complicated diagrams formed by decorating 
Fig. 2 by the diagrams of Figs. 3 and 4. Thus, the equations in the text 
describing Fig. 2 have a wider applicability. However incorporation 
of Fig. 4f would require us to replace Aj by a memory kernel in the 
equations describing the diagram. See text. 



dfj/dt = l^frUrj - l^fjUj, 



df'kldt = tf'sa'.k - ^f\a'ks 



(18a) 



(18b) 



where^ is the fraction of systems in state j and «,> is the 
rate of jumping from state j to r. Multiplying the first 
equation by f\ and the second by fj we obtain 



It can also be shown using the methods developed 
previously [38] that Fig. 4b is equivalent to Fig. 4a. 
Specifically, one can choose rate constants for the up- 
ward and downward steps in Fig. 4a that are com- 
pounded from those of Fig. 4b in such a way that the 
occupation of the bottom well in Fig. 4a equals the sum 
of those in Fig. 4b in both the equilibrium and the flux 
determined [38] steady state solutions. 



d(/"'(i/;)/df = ^fkfrOirj - '^fkfiOijr + %fjf\a\k " l^fjfkOi'ks 



S l,frf\(a,j8,k + a\k8rj) - S ^fjf\{ajAk + a'faS^,) 



SS (/;/'. )A„;,.,-SS(;5AH 



^jk;rs ■ 



(19) 
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Or, if we relabel the indices so that € = (r,s) and ; 
(j,k) then we obtain 



dN/dt = d\DN)ldx'^ - d(vN)/dx -bN + AN (22a) 



dN./dt = XNe An - XNAie 



(20) 



which is the master equation for the composite system. 
Notice that the complexion, N, , of the composite system 
is the product of the complexions, fj, of the individual 
systems, but the composite transition coefficients are 
sums of the individual transition coefficients. These re- 
sults are readily generalized by the process of induction 
to a system consisting of any number of subsystems, the 
only condition being that the subsystems do not interact 
with each other. We again see that the complexion of the 
composite system is the product of the complexions of 
the individual systems, but the composite transition 
coefficients are sums of the individual transition coeffi- 
cients. 

Thus, if we could find a smallest set of independently 
interacting molecules we will have simplified our prob- 
lem considerably. Fortunately there is a confluence of 
intuition and experiment that suggest that this can be 
done. First, what is happening at point a cannot be 
influenced by what is happening at point b provided that 
the two points are sufficiently far apart. So, there is a 
smallest size. Second, this size seems to be very small 
indeed. Stillinger, on the basis of computer modeling 
and other considerations has concluded [39] that the 
number of molecules involved in the basic diffusion step 
is on the order of several molecules for simple van der 
Waals systems. Perhaps a local density decrease allows 
a molecule to jump out of a cage, or perhaps two 
molecules interchange, resulting in a net flow. 

As a result of these considerations we can maintain 
that the Nt, Mi of Eqs. (17) refer not only to configura- 
tion space, but also to particles or quasiparticles in our 
3-d space. A connection is thus made between the 
trapping model of Di Marzio and Sanchez [32] who 
trapped the configuration point and the trapping model 
of Odagaki et al. [40] who trapped atoms. Of course 
trapping atoms implies trapping the configuration point 
and conversely. The context of the discussion easily de- 
termines what kind of particle or quasiparticle is being 
trapped. 

Equations (17) can be transformed into a continuum 
version by using 

Nj{t) -> N{t,x), Nj+, (0 -> N(t,x + Ax), Mj -> M(t,x), 
etc. 

aj -> a(x), fij -> ;8(x), Uj+i -> a{x + Ax), etc. (21) 
we obtain 



dM/dt = + bN-AM, 



(22b) 



where 



D = (Axf(a + l3)/2 and v = (Ax)(q; -;S). (23) 

D, V , b, and A can all be position dependent. 

The rate constants are determined as follows. From 
Eq. (15) we have 



Aj = bj e\pi-\Ej\/kT), bj = b. 



(24) 



where Ej is the depth of the well. We argue that the 
energy appears only as a barrier restricting the escape 
from the wells-there is no attraction of the phase point 
into a well. The bj are also all chosen to be equal be- 
cause we can think of nothing that distinguishes them 
from each other. Allowing the aj to be different from 
the Pj accounts for a drifting of the phase point towards 
a region of phase space. This should be useful if we 
impose an external field. If we assume no x dependence 
for a and 13 then D and v are constants and the d^DN/dx^ 
term is the ordinary diffusion term. Our Equations now 
read 

dN(t,x)/dt = Dd^N(t,x)/dx^ 

- vdN(t,x)/dx - bN{t,x) + A{x)N{t,x) (25a) 



dM(t,x)/dt = bN(t,x) - A(x)M(t,x), 



(25b) 



where we have written all t, x dependencies explicitly. 
Since a and 13 are much greater than b we know that 
after jumping out of a low lying well the phase point 
will travel extensively horizontally before being cap- 
tured by a deep well. Since b does not depend on x and 
is not a function of well depth the rate of filling the 
wells is random. Thus the horizontal distribution of well 
depths which we assume to be random along the chain 
(see Fig. 2) is unimportant. If W(E) is the number of 
wells of depth E then they are filled with a rate propor- 
tional to W(E). Over a large period of time the escaping 
from wells is determined by both W(E) and the rate of 
escape (exp(— ;S l^l)) from individual wells. This 
allows us to replace the distribution of wells by wells of 
one depth. In this case Fig. 2 becomes simplified even 
further so that the vertical lines have the same length. 
The equations now can be simply solved since A now 
has no x dependence. Using the method of moments on 
Eqs. (25) we find 
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d{N)IAt = -b {N)+A{M) 
d{M)IAt = b {N)-A{M) (26a) 

d{xN)ldt = v {N)-b (xN) +A{xN) 
d {xM )/dt = b {xN)-A {xM ) (26b) 

d{x'-N)/dt = 2D (N) + 2v (xN) 
-b {x^N)+A{x^M) 
d{x^M)/dt = b (x^N) -A (x'M) . (26c) 

The nice thing about these equations is that we can 
solve the nth pair of equations for the nth order mo- 
ments in terms of the lower order sets. We will exploit 
this fact in the next section. 

Finally, considering only the sequence in time of the 
occupation of the deep wells by the configuration point, 
with bj = b, ItNj = nN, and assuming that horizontal 
motion is so fast that Nj = N, the sum over the Nj in 
Eq. (17a) yields. 



ndN/dt = -nbN +2A,M, 
AM^Idt = bN - A,M, 



(27a) 



(27b) 



Here the total number of shallow wells is n. Figure 5 
displays the diagram associated with these equations. 

One notes that Eqs. (17), (22), and (25) are very 
similar to equations arising in modeling chromatogra- 
phy [41]. In that case the diffusion and drift terms 
model the behavior of the eluting material as it travels 
along in the mobile phase. Nit, x) being the amount of 
material in the mobile phase, while M{t, x) is the 
amount of material adsorbed on the adjacent surface or 
in pores [42]. 

Our minimal models are all now well defined and 
deriving their implications is merely a matter of mathe- 
matics, albeit sometimes very difficult mathematics. 
The remaining conceptual problem, to which we now 
turn, is to relate the solution of these minimal models to 
the frequency and temperture dependent complex vis- 
cosity T]*{co,T), diffusion coefficient D{co,T) and 
dielectric response e{co,T). 

2.4 Insights From Our Minimal Models: Deriva- 
tion of £> (0, J) And rj(0, r) 

2.4.1 The Diffusion Coefficient 0(0,7) When 
All Wells Have the Same Depth Equations 26 are 
easily solved for the moments. After some labor, with 



obvious assumptions on the initial conditions we obtain 
to first order in the drift velocity 

(x) = (xiN + M)) /{{N + M)) = (A/ib + A))vt (28) 

( (X - (x )f) = ( (X - (x )nN + M))/( (N + M)) 

= (A/ib + A))2Dt . (29) 




Fig. 5. If the a and /3 of Fig. 2 are very much larger than the b and 
A then we can argue that the configuration point nmning about in the 
"configurational sea" sees an unbiased statistical sample of the wells 
before falling out of the "configurational sea" into any one of them. 
Thus if we are interested only in the sequence in time of occupation 
of the wells by the configuration point the diagram of Fig. 5 suffices. 
In the text the simplified equations describing Fig. 5 are obtained. 



Notice that the diffusion coefficient is diminished by 
the factor A/{b+A) (because from Eq. (24), 
A/b = exp(- \Ej\lkT) and the wells are deep we will 
ignore the A in the denominator of A/(fo + A)). These 
equations have the obvious interpretation that every- 
thing, both drift and diffusion, is being slowed down by 
the factor A /fo which is the ratio of jump rates. As long 
as the particle is in a deep well there is no activity. Any 
resulting activity is proportional to the rate of escape, 
exp(- \Ej\lkT), from the deep wells. 

We now seek to further interpret this result. 
The ordinary diffusion equation without sinks (dN/dt = 
Dd^N/dx^) has as its Green's function the Gaussian 
distribution (477D/)"^exp(-x^/4D/). In the probabihstic 
formulation of the diffusion equation this Green's 
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function has the physical interpretation of representing a 
random walk as in Fig. 6a. There is no pausing between 
steps of the random walk. However, the equations of our 
minimal models have the interpretation that when the 
particle is in a deep well there is no motion until, after 
a long time the particle escapes the well. Thus, in the 
probabilistic interpretation of our minimal models our 
physical process is represented by a random walk with 
a pausing time between steps. The steps themselves 
correspond to the horizontal motion characterized by 
the diffusion constant D while the pausing corresponds 
to the time spent in the deep wells. Thus, the effective 
diffusion coefficient is 




aN(t,x)/at = Da^N/ax^ 

-bN(t,x) +A(x)N(t,x) 
aM(t,x)/at = bN(t,x) -A(x)M(t,x) 




Fig. 6. a) In the probabilistic interpretation of the diffusion equation 
dC/dt = DV^C the Green's function represents a random walk with 
no pausing time between steps of the random walk, b) In the 
probabilistic interpretation of the equations describing our minimal 
model the trapping in deep wells corresponds to a pausing time 
between steps of the random walk. The spatial aspects of the walks are 
identical in both cases. 



Drff(0, T) = (Ax)V2(A0eff 
= (A/b)D = ((A/b)(Axf/2At 
= exp(-| £, \/kT){Axf/2At . (30) 



2.4.2 The Viscosity rj(0, T) When All Wells Are 
of the Same Depth In Eq. (30) we have taken the view 
that the paths traversed in configuration space are the 
same for both the case of pure diffusion and that of 
diffusion with traps (See Fig. 6). This means that the 
only difference between the two cases is the time to take 
each step. For diffusion with traps we write 



Afef 



At + Afweii, At « Af„ 



(31) 



where Afweii is the time spent in the traps between 
jumps, while At is the time spent traversing the path in 
the configurational sea (the time spent between jumping 
out of one well and falling into the next well). Since 
viscosity is inverse to diffusion we will assume that the 
viscosity is proportional to the average time spent in the 
deep wells. This notion is verified in Appendix B. Thus, 

y}{Q,T)IB ^{t)= I tP(t,T)dt (32) 

•'0 

where P{t,T) is the normalized probability density that 
the configuration point escapes the well at time t. 

When all the wells are of the same depth P{t,T) is 
easily calculated. The probability ^{t,T) of the particle 
being in the well at time t is 

^{t,T) = exp{-btexp{\E\/kT)) 
= exp(- f/r); r = fe-'exp(+ \E \lkT) (33) 

and the probability density P{t,T) for exiting the well 
at time t is 



P{t,T) = d-^ldt = T-'exp(- tiT) 



(34) 



The exponential approximation for ^ is a good one. 
To see why consider the configuration point in a well. 
It decays exponentially initially. This can be seen by 
solving the generic matrix equation to which the set 27 
belongs 



dJV/df = AN 



(35) 



where A^ is the set {Nj,Mj) and A is the matrix of Eqs. 
(27). We obtain [43] 
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N=exp(+At)N(0) 



(36) 



and if we begin with one particle in one well we see that 
for small t we have our exponential decay. But when the 
particle jumps out of this well the chance that it comes 
back into the same well is very small since there are so 
many other wells. Thus, we are confident of our as- 
sumed form [(Eq. 34)]. However, it is stressed that Eqs. 
(27) should be solved rigorously to bolster the argument. 

Equation (34) when substituted into Eq. (32) gives 

r](0,T)/B '^{t)= tP(t,T)dt 
Jo 

= tT-' exp(- t/T)dt = T (37) 

•'0 

which was to be expected. 

2.4.3 ij(0,r) When the Wells Are of Different 
Depths However, solving the problem where the deep 
wells are all the same depth is not the same as solving 
the problem for glasses since glasses have a distribution 
of well depths. We need to evaluate P(f, T) for this latter 
case and also calculate a new effective diffusion coeffi- 
cient. P{t, T) is exactly calculable from Eqs. (27) since 
in the probabilistic interpretation the configuration point 
jumps from well to well, and there is no Ax involved in 
Eqs. (27). A configuration point in a well of depth E 
sees only the barrier and therefore the probability that it 
be in the well at time t is given by Eq. (33). Let W(E) 
be the weight distribution for wells of depth E. Notice 
from Eq. (37) that JW(E)P(t, T)dt = Z?-'W(£)exp(l£,/ 
kT) which states that the time spent in wells of level E 
is given by the Boltzmann factor weighted by the degen- 
eracy factor W(E). This is in perfect accord with the 
ergodic theorem. An estimate of the relaxation function 
p(t,T) describing the exiting from wells can now be 
made by weighting the distribution function P(t,T) (see 
Eq. (34)) for the occupation of the well of depth E by the 
weighting function W(E). 

pit, T) = SWiE)Pit, T)dE/SWiE)dE, 

= JW\E)P(t,E)dE (38) 

W'(E) = W(E)/JW(E)dE (39) 

The viscosity becomes 

V(0,T) <^ (t) = b-' /W(£)exp(+ \E\/kT)dE 

= JW\E)T(E)dE (40a) 



The right-hand-side of Eq. (40a) is closely related to 
the partition function. We develop the consequences of 
this in Sec. 2.6. 

In Sees. 2.4.2 and 2.4.3 we presumed that the process 
of flow could occur if only one particle jumped out of 
its well. But suppose it is required that within a space of 
a given volume there needs to be M particles that have 
simultaneously jumped out of their wells in order to 
have flow. It is shown in Appendix C that Eq. (40a) is 
generalized to 



r](0,T) a {t)" = [b-' /W(£)exp(+ \E\lkT)dEf 



UW\E)T{E)dEf 



(40b) 



This allows us to express the temperature dependence 
of 17 as 

logT](0, T) = B + Mlog[/ W(£)exp(+ \E \/kT)dE] 

(40c) 

where B and M are considered to be constants. 

2.4.4 D(0,T) When the Wells Are of Different 
Depths We now seek to calculate the diffusion coeffi- 
cient when we have a distribution in well depths. The 
answer to this can be obtained by solving Eqs. (27) or 
(17), but we are unable to do this presently. Instead, we 
argue that the diagram of Fig. 2 which is our model for 
real glasses can be approximated under certain circum- 
stances by the simpler diagram with all wells being of 
equal depth provided we choose an effective well depth. 
We choose for this effective well an effective rate con- 
stant Aeff given by 



(SWOM,ff = SW,M, 



(41) 



The form of Eq. (41) reduces to the proper limiting 
form when there is only one well depth and additionally 
allows the escape from very deep wells to be the rate 
determining steps. The Wi appear as shown because the 
number of times a particle falls into a well of depth £, 
is given by W, . The argument for this is that as soon as 
a configuration point escapes its well, because of the 
large value of D while running about in the upper wells 
it has exposed itself to the other wells, and because b is 
independent of x it falls into each well with equal 
probability. If the number of wells of depth j is W, the 
configuration point falls into a well of energy Ej with a 
probability W, and then tries to escape with a probabil- 
ity proportional to Aj. Thus, we know that Wj is propor- 
tional to the number of well of type j and the effective 
diffusion coefficient for Fig. 2 is then given by 



D,, 



■■ AafDIb 



(42) 
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2.5 Evaluation of the Frequency Dependent 
rj(fe>,r), e{o),T)andD{o),T) 

2.5.1 Evaluation of r}(o},T) Equation (10b) has 
its analogue in polymer physics The complex viscosity 
is 

r]*{w,T) = G{w,T) = exp{-iwt)g{t,T)dt (43) 

•'0 

and the frequency dependent shear modulus is defined 



G*{co,T) = icori*{co,T) . 
At zero frequency we showed that 



(44) 



-i 



r]*{0,T) ^{t)= I {tT-')exp{-t/T)dt = t . (45) 

•'0 

But it would be wrong to identify g{t,T) with the inte- 
grand of Eq. (45). In fact since 



I 



(lM!)(fT-')"exp(-f/T)df= T 



(46) 



any value of n would be permitted if the sole criterion 
were that the integral equal t. Formulated in this way it 
is obvious that n = gives the correct g{t,T) since it 
corresponds to a Maxwell element. Thus g{t,T) is pro- 



1 h 
Jo 



portional to | tP{t, T)dt and since the value of g (0, T) is 
Go we have 

g{t,T) = GoSxp{-t/r). (47) 

This gives immediately 

r]*{w,T) = G{w,T) = Got/(1 + iwr) (48) 

17* = rj' - irj" , (49) 

T?' = Got/(1 + w\^) (50a) 

T?" = Go«tV(1 + «'t') (50b) 

while for a distribution W'{E) of well depths we obtain 

r]*{w,T) = JW\E)Got/{1 + iwT)dE (51) 

g{t,T) = JW'{E)GoCxp{-t/T)dE (52) 

Ty' = JW\E)Go t(1 + «V)d£ (53a) 

T?" = JW\Ey{l + w\^)dE . (53b) 



These relationships show clearly that non-Debye 
frequency behavior occurs because there is a distribution 
of relaxation times. 

2.5.2 Evaluation of Dielectric Response 
e{oi, T) Granted the calculation of the complex vis- 
cosity, the dielectric constant s{co,T) can also be ob- 
tained. Debye showed that if the dipoles are each imag- 
ined to be imbedded in the center of spheres (one dipole 
per sphere) that are in turn imbedded in a viscous fluid 
of viscosity rj then the dielectric response is easily cal- 
culated [44]. Based on this result Di Marzio and Bishop 
showed [45] that if the viscous fluid has a complex 
viscosity ri*{co,T) then the formula is a simple general- 
ization of the Debye formula, the only change being that 
ri*{co,T) replaces 17(0, T). Thus, 



[s(c,T) -s(cc,T)]/[s(0,T)-s(cc,T)] 



(1 + iM'r]*(M,T)Ay 



(54) 



where A is a dimensional constant. The plus sign occurs 
in Eq. (54) because of our choice of the convention for 
the Fourier transform as in Eq. (43). This is consistent 
with Ferry's [46] development of viscoelasticity for 
polymers. 

2.5.3 Evaluation of D(m,T) Equation (25a) 
shows that the diffusion coefficient D is a constant. In 
order for it to have a frequency dependence we would 
have to have had jD{t-T)d^N/dxMT for the first term 
on the right hand side of Eq. (25a). But, this is not the 
case. Equivalently we could have used the folding opera- 
tion and written D{t-T) = Dh{t- t). Further, from Eq. 
(29b) we see that the effective diffusion coefficient also 
has no frequency dependence at least to the quadratic 
approximation. Therefore, for our model we expect no 
frequency dependence in the diffusion coefficient. 



D{co,T) = b-' exp{-l3\E\)D 



(55) 



For a distribution of well depths we have as before 

Dico,T) = D JWiE)dE/[JWiE)expi+ l3\E\)dE] . 

(56) 

2.6 Evaluation of W{E) 

The above relationships are quite remarkable for they 
state that long time relaxations-viscosity, diffusion and 
dielectric response depend only on the well depths and 
the distribution of well depths. The only thing remain- 
ing is for us to evaluate W{E). Notice that if this can be 
done then our kinetics of glasses will depend only on the 
equilibrium statistical mechanics. For glasses statistical 
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mechanics plus the principle of detailed balance is 
everything provided we are looking only at the long time 
behavior. 

The classical and quantum mechanical partition func- 
tions are given by (we ignore the thermal wavelength) 

Gciassicai = JsKp(-l3E(. .qj . .))TTdqj (57a) 

Gqm. = Jexp(-i3£,)d, (57b) 

where the integral signs represent discrete sums and/or 
continuum integrals. By grouping together all states 
with the same energy we obtain 

Gciassicai = J sxp(- l3E)W(E)dE (58a) 

Gqm. = Jsxp(-l3E)W(E)dE (58b) 

which are identical in form to the argument of the 
logarithm on the RHS of Eq. (40c). Using the formula 
f c = -A:nnG which connects the configurational part of 
the Helmholtz free energy F^ to the partition function Q 
we have immediately 



logT?(0,r) = B - M'FJkT . 



(59) 



This remarkable formula which relates viscosity to free 
energy is very different from the Vogel-Fulcher- 
Tammann-Hesse form [47], the Bendler-Shlesinger 
form [48], the Avramov form [49], the Adam-Gibbs 
form [12] or the mode coupling theory result [50]. We 
discuss it in Sec. 3.1. 

The frequency dependent viscosity, given by 
Eqs. (53), cannot be expressed as a function of free 
energy. Rather, we first must determine W{E) sepa- 
rately before we can evaluate ri(w,T). If in Eq. (58a) we 
choose the lowest energy as our zero of energy, then 
exp(- I3F(I3)) is the Laplace transform of W(E) and 
W(E) is the inverse transform of exp(-^f (^)). 

Another approach is to use the results of Stillinger 
who suggests that W(E) is given by [51, 52] 



W(E) = exp(- 0(E-Eof). 



(60) 



With this substitution the time dependent shear mod- 
ulus, Eq. (52), reads 

git,T) = Sexpi-eiE-Eof) 

exp(-btexp(-l3\E\)dE I Jexp(-0 (£'-£'o)')d£' . (61) 

The time dependent behavior of Eq. (61) is closely 
related to that of the "after-effect function" tabulated by 
Janke and Emde [53]. As shown previously the after- 



effect function has a time dependence which looks very 
much like the stretched exponential function [23]. In 
fact, Stillinger, starting from the empirically observed 
stretched exponential form for relaxation shows that the 
Gaussian form for W{E) is implied [52]. 

3. Discussion and Conclusions 
3.1 Discussion of Results 

Equation (59) which connects viscosity to free energy 
is remarkable in several respects. First, it states that the 
viscosity and its temperature derivative are continuous 
as we proceed through the transition. We had, in Sec. 
1.8, used the argument that the dissipative quantities 
should have the same transition behavior as the thermo- 
dynamic variables. So, for a first-order transition the 
viscosity is discontinuous through the transition because 
the entropy and volume are. But we have now obtained 
the result that for a second-order transition the viscosity 
does not show a break as we traverse the transition point. 
In the past various groups have argued that the volume 
[54] is the controlling quantity, or the enthalpy [55], or 
the entropy [1-5]. We are claiming that the entropy 
theory of glass formation, which is merely a theory that 
locates the transition in temperature and pressure space 
as a function of the molecular parameters such as chain 
length, intermolecular energies and intramolecular stiff- 
ness energies etc. (see Sec. 1.4) can be extended to 
include slow motion kinetics. When this is done the only 
determinate of the kinetic aspects of glass formation in 
the limit of zero frequency is the thermodynamic free 
energy! See Eq. 62e. However, as Eqs. (51-54) show 
this is not true for the frequency dependent dissipative 
quantities. 

The Vogel-Fulcher-Tammann-Hesse form [47] from 
which the WLF equation [56] is easily derived is 

logT? = B+A/(r-ro). (62a) 

The Bendler-Shlesinger form [48] is 

logT? = B + A/(r-ro)'' (62b) 

The Avramov form [49] is 

logT? = B + 0A3A{AITy (62c) 

and the Adam-Gibbs form [12] is 

logT? = B +A/TS, . (62d) 

These forms should be compared to our form which 
is 

logT? = B- AFJkT. (62e) 
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We will not discuss the mode coupling form for viscos- 
ity since we accept the argument [57] that the implied 
singularity is considerably higher than T^. 

Although each of the first four forms has some theo- 
retical underpinning it is probably true that the reason 
they fit experimental data well is that they (the first 
three) are three-parameter fits and the viscosity curves 
are rather structureless to begin with. To see that it is not 
terribly significant to fit a curve of relatively little struc- 
ture with three parameters imagine B to locate the curve 
vertically, another of the parameters stretches the curve 
so that there is a fit at both high and low temperatures. 
Finally the third parameter gives the curve the proper 
amount of curvature. Viewed in this way we see that the 
fact that formulas of different construction give decent 
fits to the data is not surprising. A real test of the 
theories is whether they can determine the values of the 
three parameters from theory. 

Viewed from this perspective the last two Eqs. (62d), 
and (62e) are more significant because they contain one 
less parameter. The original GD lattice theory can be 
used to obtain F^. A real theory should contain no 
parameters. Schroedinger's equation plus the laws of 
statistical mechanics should be sufficient. The authors 
intend to examine the meaning of the B and A parame- 
ters of Eq. (62e) in a subsequent paper. For now we will 
merely comment on the implication of the form of our 
equation, assuming A and B to be temperature indepen- 
dent. 

Angell's classification [7] of glasses into strong and 
fragile receives an easy interpretation from Eq. (62e). 
First, we need to use the experimental value of the free 
energy in Eq. (62e). There is a general consensus that 
the specific heat break at the glass transition, Cp,c, varies 
inversely with temperature [58]. We therefore use the 
form Cpc = a/T. 

Cp,c = a/T -^S,= a{l/T2 - \IT) -^ (63) 

F, = -C-a {TIT2 -\) + a \n(T/T2l T2 < T (64a) 



-C, r < r. 



(64b) 



where the constant of integration C is (part of) the 
energy of activation. 

To obtain these equations we integrated Cp^ = TdSJ 
dT, S^ = dFJdT and ignored any pressure dependence. 
Below The transition temperature T2 the configurational 
entropy is zero according to the simple version of the 
GD theory so that we have only energy of activation 
while above T2 the specific heat is assumed to decrease 
inversely with temperature in accord with experiment. 

Using Eq. (62e) we can eliminate B by choosing a 
reference temperature T* for which the viscosity equals 



10'^ poise. A little algebra results in 
logT] =13 + ^xlnx + (l-x) 
X (^[1 +ln{T*/T2)]-0) 

dlogrj/dx = 6-^\n{T*IT2) + ^\rvc 



(65a) 



(65b) 



a'logT]/ax' = ^Ix, T2<T*< T, 

T2<T<T* (65c) 

logT] = 13 + e{x-\)-^ [{TVT2 - 1) 

+ ln{T*/T2)] (66a) 

logrj/dx = 9 (66b) 

dHogr}/dx2 = 0,T<T2<T* (66c) 

logT] = 13 + e(x - 1) + ^[(T*/T2 - x) 

+ xlnixTJT*)] (67a) 

dlog-q/dx =6-^ \n(T*/T2) + ^\wc (67b) 

d^logrj/dx^ = ^/x, T* <T2<T (67c) 

logT] =13 + e{x - 1) (68a) 

dlog-^ldx = e (68b) 

d^\ogT]ldx'^ = 0, T ^ T* <: T2, 

T* <T< T2 (68c) 

where 6 = CA/kT*, ^ = oA/kT*, x = T*/T. T* is the 
temperature for which 17 = 10" poise. If we had picked 
10^ as the reference viscosity then the above equations 
would be the same with y replacing 1 3 and T* being the 
temperature at which the viscosity is 10^ poise. 

Equations (65c,) and (67c) show that the curvature is 
positive (curve is concave up) and that the curvature is 
greater the larger the specific heat. Also, as the value of 
T*/T decreases the curvature is larger. Below the glass 
temperature we predict pure Ahrennius behavior. These 
features are also features of Angell's classification of 
glasses into strong and fragile varieties. An interesting 
prediction is that if T*/T2 = 1 then the initial slope at 
T*/T= 1 is independent of specific heat. It does however 
depend on C. 
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We can test these predictions for polymers using data 
for polydimethylsiloxane of varying molecular weight. 
Roland and Ngai [59] using dielectric relaxation data of 
Kirst et al. [60] and specific heat data of Bershstein and 
Egorov [61] created fragility plots of the logarithm of 
relaxation time versus TJT where T^ was defined as the 
temperature for which the relaxation time was one sec- 
ond. These curves which are reproduced in Fig. 7 show, 
as Roland and Ngai observed, 1) that the slope of the 
curves at T*/T= 1 are independent of specific heat — we 
predict this, 2) The curvature is larger the smaller the 
value of Tg/r — we predict this, and 3) the curves flare 
out for low Tg/r with the higher specific heat (low 
molecular weight) material flaring up and the low 
specific heat (high molecular weight) material flaring 
down — we predict this. The filled circles are our numer- 



ical predictions. We chose A and B to fit the center 
curve. We then scaled ^by the ratio of the specific heats 
for the low and high molecular weight polymers to ob- 
tain the upper and lower points at each temperature. Our 
fits assume that C is independent of molecular weight. 
We also give the formulas for the case that the config- 
urational specific heat is constant above T^. Our reason 
for doing this is that although the GD lattice model 
predicts that the configurational specific heat ap- 
proaches zero as the temperature increases it does not 
do so with purely inverse temperature dependence. So, 
a combination of the two specific heat variations may 
better fit the experimental data. 
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Fig. 7. An attempt to explain the fragility plots of Angell. Using for the configurational free 
energy a form derived by assuming that the specific heat is proportional to T"', which is in 
accord with experiment, we obtain a fit to the plots of log relaxation time versus TJT ■ The 
curves are experimental data for polydimethylsiloxane of varying molecular weight, and the 
circles are predicted values. That 1) the curves all start with the same slope at T*IT = 1; 2) the 
curvature increases with decreasing T*IT; 3) the curvature increases with increasing specific 
heat are all predicted by our equation. See text. 
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S, = a'ln(T/T2), Tj < T (70a) 

S, = 0, r < 7^2 (70b) 

F, = - a'Tln(T/T2) + a\T-T2) - C, 7^2 < T (70a) 

F, = -C, T<T2 (70b) 

If we again define T* as the temperature for which 
the viscosity equals 10'^ poise we obtain 

logT? = 13 -rin(x) + {CT2IT* + e){x - 1) (71a) 

d\ogT]ldx = CT2IT* + B-Cx-\ (72b) 



dHogrj/dx^ = + ^'x ^, 

T'2 s r* < r T'2 s r < r* 



(72c) 



logT? = 13 - ^'ln(r*/r2) + Bix - 1) 

+ ^'(1 - T2IT*), (73a) 

d\ogT]ldx = e, (73b) 

dHog-q/dx^ = 0, r < 7^2 < r* (73c) 

logT? = 13 + ^'in(r/r2) + eix - i) 

+ f(T2/T*)(x - T*/T2) (74a) 

dlog-q/dx = C'TiT* + e- a'x-', (74b) 

dHogT]/dx^ = + ^'x^ , r* < 7^2 < r (74c) 

log'r?= 13 + 0(x- 1), (75a) 

logTj/to = e, (75b) 

d^logrj/dx^ = 0,T* <T <T2,T <T* <T2 

(75c) 

where ^= a'Alk and B = CA/kT*. These curves again 
show the features of the strong-fragility plots discussed 
by Angell. 

It should be noted that if either of the above forms for 
the entropy is substituted into the Adam-Gibbs form 
[Eq. (62d)] one obtains a decreasing slope with 
increasing specific heat at T*/T= 1. Also the curvature 
of the logi7 vs T*/T curve becomes smaller as T*/T 



decreases which is contrary to the sense of virtually all 
experimental results. 

Can T* ever be less than T'2? Under the paradigm of 
the Vogel-Fulcher equation this is a foolish question. 
However, since the viscosity and its derivative are, 
according to Eq. (62a), continuous through the second- 
order transition and since the viscosity is never infinite 
T2 can not be located accurately by measurements of 
viscosity; we see no reason why it can not be greater 
than T*. The possibility that T2 corresponds to a finite 
viscosity may well be masked by the process of falling 
out of equilibrium which can be discussed only by 
examining the time or frequency dependent viscosity. 

The new formulas for viscosity suggests several new 
directions. First, an examination of the way which C/T* 
varies with material should be made. 

We remark that these questions are equilibrium ther- 
modynamic and statistical mechanical questions so that 
their investigation should not be difficult. For systems 
with constant C the initial slope of the curve at T*/T = 1 
would be inverse to T*. Also systems for which the 
motion is highly cooperative would show a higher C. 
Systems which have the same scaled potential energy 
surface, i.e., hE(. . qi . .) where h is any constant, should 
display superposed fragility plots. Such systems which 
have no specific heat break at T^ should all superpose 
with the form of a straight line. Finally how the specific 
heat a relates to CIT* should be examined. 

Another possibility that deserves serious consider- 
ation is that the parameters B and A have a temperature 
dependence which must be added to that of the free 
energy. This thought is consistent with the view ex- 
pressed by some that the temperature dependence of 
viscosity and diffusion at higher temperatures is ade- 
quately addressed by mode coupling theory and that the 
behavior over the full temperature range can be obtained 
by a cross-over treatment that combines the high tem- 
perature mode coupling theory with a theory of low 
temperatures such as has been presented here. 

We leave such a development to the future. 

3.2 Conclusions 

This paragraph describes the logic of our develop- 
ment. We first observed that there must exist at low 
temperatures an equilibrium glass phase because the 
crystal phase is not ubiquitous. It is only for systems that 
can crystallize that the glass phase can be considered 
to be a metastable phase. We next showed that the 
Gibbs-Di Marzio (GD) theory [1-5] which postulates 
that the glass transition occurs when the configurational 
entropy approaches zero locates the glass transition cor- 
rectly in temperature-pressure space for a wide variety 
of experiments. It also resolves the Kauzmann paradox 
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[6]. We next observed that the correct equihbrium the- 
ory of those materials that form glasses, whatever it may 
be, must be used as a groundform onto which a proper 
kinetic theory of glasses is constructed. The connection 
between kinetics and equilibrium was then made via the 
principle of detailed balance which relates the ratio of 
the rates for jumping to and from a pair of states to the 
free energy difference between the states. This law when 
combined with the observation that the configuration 
point of a glass system spends most of its time in deep 
potential energy minima allows us to construct a mini- 
mal model (a trapping model) which can be solved in 
some limiting cases. We thereby obtained formulas for 
the complex viscosity ri*(co,T) and shear relaxation 
modulus g(t,T), the diffusion coefficient D(co,T) = 
D(0, T) and the dielectric response e(co, T). 

Our relaxation modulus has the form of the after-ef- 
fect function tabulated by Jahnke and Emde [53] which 
we had obtained previously [32]. Its behavior is very 
close to the stretched exponential form. 

More surprising is our formula relating the zero fre- 
quency viscosity 17 to the configurational part of the 
thermodynamic Helmholtz [62] free energy F^ 

log 77 = B - AFJkT 

It is surprising that the viscosity is continuous 
through the transition. If this conclusion holds, viscosity 
or other dynamical measurements may be the worst way 
to locate glass temperatures. The use of thermodynamic 
quantities which show breaks in slope should be pre- 
ferred. Initially we had expected (naively in retrospect) 
that there should be a break in slope of 17(7") vs T. 
Below the transition the behavior is Arrhenius corre- 
sponding to the fact that there is energy but not entropy 
of activation below the transition. Above the glass tran- 
sition entropy of activation kicks in. 

When sensible approximations for F^ are used this 
formula displays the main features of the strong-fragile 
glass classification scheme proposed by Angell [7]. 
Glasses with small specific heat breaks at the glass 
transition show little curwture on log 17 versus inverse 
temperature plots while glasses with large specific heat 
breaks show positive curvature. See Sec. 3.1. 

We have not yet examined the temperature depen- 
dence of B and A in the above equation. An approach to 
this problem is to excise from phase space those phase 
points corresponding to deep potential energy minima 
and solve the kinetics of such a circumscribed space. 
Since the resulting equations should be applicable to the 
high temperature side of the glass transition it may be 
that mode-coupling theory can be used for this part of 
the problem. A theory of glasses that is valid over a wide 



range of temperatures undoubtedly requires incorpora- 
tion of vibrational properties. 

4. Appendix A: A Particle Falling Back 
into a Well From Which It Is Trying to 
Escape Retains Its Exponential Dis- 
tribution 

Let a particle have a normalized probability distribu- 
tion p(t) for escaping from a well. But after it has left 
the well let there be a probability 1-0 that it fall back 
into the well and cfi that it escape permanently on that 
attempt. Then the particle can escape permanently after 
it left for the first time, or the second time or the third 
time, etc. The true distribution function for escape is 



Ptm. = p(t)(j) + S pitl)i\-4>)pit2)dttdt2(t) 



+ S Piu)p{t2){\-4)fp{h)dhdhdh4) +■■■ 
ti+t2 + h = t 

(Al) 

The integrals are recognized to be folding operations. 
Denoting the Laplace transform of p by P and taking 
the Laplace transform of Eq. (Al) we obtain 

+ p\i-(t)f4> + . . . 

= P0/(1_(1_0)P) (A2) 

This equation states that if p is normalized, so that 
^(0) = j p(t)dt = 1), then /^tme is also normalized. 

For the special case where p is exponential in t we 
have 

pit) = T-' exp(-f/T), P = Pis) = T-'/s + T-') . (A3) 
Using Eq. (3) in Eq. (2) we obtain 

P^. = P^ru. is) = <i>r-'lis + <i>r-') (A4) 

which on transforming back to the time domain gives 

/"true = 0T-' eXp(- (i>tlT) (A5) 

Thus, the final distribution remains exponential and 
is normalized but the time constant for exiting the well 
is increased by the factor <^"\ 
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5. Appendix B. Calculation of the Shear 
Viscosity r/(0,r) 

To calculate the shear viscosity consider a material 
between two parallel plates one of which is fixed and 
one of which moves under a shear force 



But 



o-/(v) = 77(0,7) 



(v) = xltf = V{U - ti)/tf = VAt/tf 



(Bl) 



(B2) 



where V is the velocity that the plate has for the time 
interval tf - ti = At during which the system is flowing 
in the configurational sea. From time zero to tj there 
was no motion because the particle was in one of the 
deep wells. We obtain 

77(0, T) = (o-(0, T)/V)tf/At = T]jf/At 

= Tjcft + At)/At = 77^(1 + ti/At) . (B3) 

This process of sticking and slipping is imageined to 
happen over and over again. The distribution function 
P{ti, T) for jumping out of a well at time t; when all the 
deep wells have the same depth is given by 

P(ti,T) = T-' exp(- t/T), = b"' exp(+ I3\E\) . (B4) 

Thus over a period of time sufficiently large to allow 
for many visits to the deep wells we have 

T?(0, T) = T?,(l + (t) I At) - {r)JAt){t) 

= WAt)T (B5) 

which proves our contention that the viscosity is propor- 
tional to the average residence time spent in deep wells. 
Evidently, if there is a distribution of well depths, 
W'{E), normalized so that the / W\E)dE = 1, we obtain 

77(0, T) = {rjJAt) J J W'{E)tP{t, T)dEdt = 

{rjJAt) J J W{E)b-'exp{+ ;S|£|)d£'// W{E)dE . 

(B6) 

Given the exponential character of P{t, E) we have 
succeeded in relating the zero frequency viscosity 
77(0, T) to the two integrals which are purely equi- 
librium quantities. 



6. Appendix C. What if Flow Requires 
Several Particles to Be Out of Their 
Wells Simultaneously 

It may be unreasonable to suppose that flow can occur 
in the region between two parallel plates when one 
particle only is out of its well. Certainly, as the amount 
of the material between the plates is increased it is more 
reasonable to expect that flow requires that the number 
of particles simulaneously out of their wells be propor- 
tional to the volume of material. We propose that the 
number required per unit volume be some large number 
M. 

Let us begin by considering the case where there is 
flow only if two particles have simultaneously escaped 
the wells. Consider one particle jumping out of a well at 
time ti, cruising the configurational sea for a time inter- 
val At = tf-ti, falling into a well, and then starting the 
process all over again. One can imagine these time inter- 
vals placed stochastically on the positive infinite half 
line. Obviously flow for the system will occur only when 
there is an overlapping of the At's of one particle with 
the A? 's of the other. The fraction of time that these time 
intervals overlap is obviously given by (nAt/Xit, + At)y. 
We imagine the particles to have jumped out of the wells 
n times, where n is very large. Refering to appendix B 
we calcuate the average velocity to be (v) to be 



(v) = V(nAt/(Mi + At)f 



(CI) 



The viscosity is 



■q = o-/(v) = (a/V) ■ [&i + nAt)/nAtf 

= r]dl+Xti/nAtf . (C2) 

Now if we require M particles to be simultaneously 
out of the wells in order to have flow we need only 
realize that the probability for this is 



Thus 



[nAt/X{ti + At)f 



T?='r?c[l +XtJnAtf 



(C3) 



(C3) 



Now if M or more particles are needed to be out of the 
wells simultaneously we have for the probability (frac- 
tion of time) for this to occur 



[1 + XtJnAtf +[\+ XtJnAtf 



+ [1 +XtJnAt]''*^ + 



(C4) 
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which leads to 

77 = 'r7,[l + Xti/nAtr [1 + nAt/Xti] 

-77JI +2/i/nA/f . (C4) 

We now need to evaluate 

Xti = nj /T"'exp(-//T)d/ = «T. (C5) 

But if there is a distribution of well depths then 



„ _ SWiE)tT-' expi-t/T)dtdE 
^'' " " / W{E)dE 



8. References 



jW{E)T{E)dE , , 
" Sw[E)dE = "<^> 



SWiE)b-'expi+ l3\E\)dE 
/ W{E)dE 



(C6) 



Therefore 



->7c[l + / W(E)b-'cxp(+ l3\E\)dEf 
"^ ~ / W{E)dEAt 



VM WiE)b-'expi+ l3\E\)dEf 
J W(E)dEAt 



(C7) 



By taking the logarithm of this equation we can cast 
it into a form usually used to compare with experiments 
and those equations created to explain experiments such 
as The Vogel-Fulcher law [47], the Bendler-Shlesinger 
law [48], and the Avramov law [49]. 

\ogr] =B + MlogiS WiE)exp (+ I3\E \)dE) (C8) 



B = log-r;, - Mlogib J W(E)dEAt) 



(C9) 



In this paper we treat B as a constant in order to focus 
on the temperature dependence of the second term on 
the RHS of Eq. (C8). Discussion of the temperature 
dependence of B is reserved for future work. 
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Note added in proof: 

Our diagrammatic approach allows us in Eqs. (30) or 
(42) to express the effective diffusion constant as the 
product of a term for jumping out of the wells and the 
diffusion coefficient for roaming about in the upper 
wells (see Fig. 2). But suppose now that the upper wells 
have the same kind of structure and that we can repre- 
sent D itself as a product of a term for jumping out of 
wells times a diffusion coefficient for roaming in the 
upper wells of the new diagram which has the same 
structure as Fig. 2. Done repeatedly we obtain a nesting 
of diagrams within diagrams and the result for the effec- 
tive diffusion constant is a product of exponential terms 
times the final D . Similarly the effective viscosity is a 
product of a partition function times the viscosity of the 
upper wells of the diagram of Fig. 2. But this latter 
viscosity is again a product of a partition function and 
viscosity, and so on for all subsequent nested diagrams. 
This leads to a linear combination of energy U and 
entropy TS^ which is not however that unique linear 
combination that results in free energy. Thus, the resul- 
tant viscosity \ogr\ = B - JS^ + KU/kT is continuous 
across the transition because U and S^ are, but the 
derivative of viscosity is discontinuous because the vis- 
cosity is no longer given by a free energy as in Eq. 
(62c). This is consistent with the fluctuation-dissipation 
theorem mentioned in Sec. 1.8 which predicts that the 
temperature and pressure behavior of the viscosity is 
continuous across the transition but the derivatives of 
viscosity are discontinuous, since the underlying ther- 
modynamic transition is second order in the Ehrenfest 
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